TABLE 2 

12 Boo V^ Dataset, Best-Fit Model Predictions, and Residuals 



MJD I Cal Date/UT WL (/im) Meas V'^ Error Pred V'^ Resid u (m) v (m) HA (hrs) 



50985.1681 


6/21/98 04:01 


2.20 


0.300 


0.022 


0.296 


0.004 


-36.90 


-101.78 


0.01 


50985.1847 


6/21/98 04:25 


2.20 


0.436 


0.027 


0.415 


0.021 


-30.47 


-103.27 


0.41 


51007.1544 


7/13/98 03:42 


2.20 


0.745 


0.039 


0.774 


-0.028 


-18.13 


-105.21 


1.13 


51007.1652 


7/13/98 03:57 


2.20 


0.750 


0.037 


0.762 


-0.012 


-13.49 


-105.67 


1.39 


51008.1611 


7/14/98 03:52 


2.20 


0.940 


0.100 


0.870 


0.070 


-14.06 


-105.62 


1.36 


51008.1747 


7/14/98 04:11 


2.20 


0.972 


0.093 


0.915 


0.057 


-8.15 


-106.02 


1.68 


51009.1677 


7/15/98 04:01 


2.19 


0.509 


0.044 


0.439 


0.070 


-10.01 


-105.92 


1.58 


51009.1871 


7/15/98 04:29 


2.19 


0.734 


0.080 


0.587 


0.148 


-1.44 


-106.22 


2.05 


51009.1888 


7/15/98 04:31 


2.19 


0.647 


0.062 


0.598 


0.048 


-0.73 


-106.22 


2.09 


51237.4622 


2/28/99 11:05 


2.20 


0.757 


0.137 


0.697 


0.060 


-42.40 


-100.17 


-0.35 


51237.4895 


2/28/99 11:44 


2.20 


0.684 


0.143 


0.698 


-0.015 


-32.20 


-102.90 


0.31 


51237.5166 


2/28/99 12:23 


2.20 


0.697 


0.102 


0.693 


0.004 


-21.12 


-104.84 


0.96 


51237.5444 


2/28/99 13:03 


2.20 


0.783 


0.055 


0.684 


0.099 


-9.15 


-105.97 


1.63 


51237.5727 


2/28/99 13:44 


2.20 


0.704 


0.063 


0.666 


0.038 


3.32 


-106.19 


2.31 


51238.3866 


3/1/99 09:16 


2.20 


0.450 


0.041 


0.442 


0.008 


-62.79 


-89.75 


-2.10 


51238.4299 


3/1/99 10:19 


2.20 


0.573 


0.041 


0.551 


0.022 


-52.03 


-96.43 


-1.06 


51238.4676 


3/1/99 11:13 


2.20 


0.767 


0.048 


0.692 


0.075 


-39.47 


-101.06 


-0.15 


51238.5350 


3/1/99 12:50 


2.20 


0.964 


0.068 


0.930 


0.034 


-12.07 


-105.78 


1.47 


51238.5531 


3/1/99 13:16 


2.20 


0.932 


0.071 


0.955 


-0.023 


-4.12 


-106.17 


1.90 


51239.3645 


3/2/99 08:44 


2.20 


0.133 


0.010 


0.125 


0.008 


-66.13 


-86.41 


-2.57 


51239.4151 


3/2/99 09:57 


2.19 


0.084 


0.010 


0.079 


0.005 


-55.44 


-94.70 


-1.35 


51239.4702 


3/2/99 11:17 


2.19 


0.133 


0.010 


0.121 


0.012 


-37.50 


-101.61 


-0.02 


51239.5062 


3/2/99 12:08 


2.19 


0.255 


0.030 


0.280 


-0.025 


-23.23 


-104.54 


0.84 


51239.5640 


3/2/99 13:32 


2.19 


0.669 


0.027 


0.701 


-0.032 


1.90 


-106.21 


2.23 


51240.5061 


3/3/99 12:08 


2.19 


0.137 


0.010 


0.149 


-0.012 


-22.10 


-104.71 


0.91 


51240.5382 


3/3/99 12:54 


2.19 


0.269 


0.010 


0.272 


-0.003 


-8.30 


-106.01 


1.68 


51240.5664 


3/3/99 13:35 


2.19 


0.419 


0.013 


0.417 


0.002 


4.17 


-106.17 


2.36 


51256.3226 


3/19/99 07:44 


2.21 


0.419 


0.025 


0.363 


0.056 


-65.43 


-87.21 


-2.46 


51256.3500 


3/19/99 08:23 


2.21 


0.446 


0.018 


0.391 


0.055 


-60.12 


-91.81 


-1.80 


51256.3980 


3/19/99 09:33 


2.21 


0.458 


0.015 


0.447 


0.011 


-46.61 


-98.71 


-0.65 


51256.4316 


3/19/99 10:21 


2.21 


0.478 


0.024 


0.494 


-0.016 


-34.54 


-102.36 


0.16 


51256.4559 


3/19/99 10:56 


2.21 


0.515 


0.018 


0.526 


-0.012 


-24.82 


-104.30 


0.75 


51256.4827 


3/19/99 11:35 


2.21 


0.556 


0.014 


0.561 


-0.005 


-13.44 


-105.67 


1.39 


51256.5048 


3/19/99 12:06 


2.21 


0.585 


0.016 


0.585 


0.000 


-3.77 


-106.18 


1.92 


51262.3474 


3/25/99 08:20 


2.21 


0.774 


0.020 


0.799 


-0.025 


-56.75 


-93.97 


-1.47 


51283.2720 


4/15/99 06:31 


2.22 


0.104 


0.010 


0.103 


0.002 


-61.07 


-91.12 


-1.90 


51283.2859 


4/15/99 06:51 


2.22 


0.095 


0.010 


0.085 


0.011 


-57.82 


-93.32 


-1.57 


51283.2972 


4/15/99 07:07 


2.22 


0.085 


0.010 


0.075 


0.010 


-54.84 


-95.03 


-1.30 


51283.3097 


4/15/99 07:25 


2.22 


0.081 


0.010 


0.073 


0.008 


-51.22 


-96.80 


-1.00 


51283.3231 


4/15/99 07:45 


2.22 


0.084 


0.010 


0.083 


0.001 


-47.00 


-98.56 


-0.67 


51283.3515 


4/15/99 08:26 


2.22 


0.169 


0.010 


0.166 


0.003 


-36.96 


-101.76 


0.01 


51283.3889 


4/15/99 09:19 


2.22 


0.412 


0.010 


0.407 


0.005 


-22.04 


-104.71 


0.91 


51316.2033 


5/18/99 04:52 


2.19 


0.067 


0.010 


0.077 


-0.010 


-55.86 


-94.47 


-1.39 


51316.2252 


5/18/99 05:24 


2.19 


0.062 


0.010 


0.073 


-0.011 


-49.50 


-97.55 


-0.86 


51316.2465 


5/18/99 05:54 


2.19 


0.086 


0.010 


0.096 


-0.010 


-42.40 


-100.17 


-0.35 


51319.2974 


5/21/99 07:08 


2.19 


0.560 


0.079 


0.660 


-0.101 


-19.17 


-105.09 


1.07 


51319.3234 


5/21/99 07:45 


2.19 


0.597 


0.020 


0.586 


0.011 


-7.91 


-106.03 


1.70 


51319.3455 


5/21/99 08:17 


2.19 


0.535 


0.030 


0.517 


0.018 


1.87 


-106.21 


2.23 



TABLE 2— Continued 



MJD I Cal Date/UT WL (/im) Mcas V^ Error Pred V^ Rcsid u (m) v (m) HA (hrs) 







51326.1502 


5/28/99 03:36 


1.62 


0.631 


0.031 


0.622 


0.008 


-62.01 


-90.39 


-2.01 


51326.1546 


5/28/99 03:42 


1.62 


0.626 


0.020 


0.608 


0.018 


-61.07 


-91.12 


-1.90 


51326.1650 


5/28/99 03:57 


1.62 


0.565 


0.018 


0.567 


-0.002 


-58.68 


-92.78 


-1.65 


51326.1751 


5/28/99 04:12 


1.62 


0.520 


0.015 


0.521 


-0.001 


-56.11 


-94.33 


-1.41 


51326.1820 


5/28/99 04:22 


1.62 


0.484 


0.016 


0.487 


-0.003 


-54.22 


-95.35 


-1.24 


51326.1890 


5/28/99 04:32 


1.62 


0.453 


0.019 


0.450 


0.003 


-52.20 


-96.35 


-1.08 


51326.2008 


5/28/99 04:49 


1.62 


0.400 


0.015 


0.383 


0.017 


-48.59 


-97.93 


-0.79 


51326.2128 


5/28/99 05:06 


1.62 


0.302 


0.012 


0.312 


-0.010 


-44.62 


-99.43 


-0.50 


51326.2192 


5/28/99 05:15 


1.62 


0.265 


0.010 


0.273 


-0.007 


-42.38 


-100.18 


-0.35 


51326.2260 


5/28/99 05:25 


1.62 


0.229 


0.010 


0.233 


-0.004 


-39.95 


-100.93 


-0.18 


51351.1581 


6/22/99 03:47 


2.20 


0.080 


0.010 


0.089 


-0.009 


-39.84 


-100.96 


-0.18 


51351.1788 


6/22/99 04:17 


2.20 


0.133 


0.010 


0.145 


-0.012 


-32.00 


-102.95 


0.32 


51351.1996 


6/22/99 04:47 


2.20 


0.247 


0.010 


0.245 


0.001 


-23.56 


-104.49 


0.82 


51351.2205 


6/22/99 05:17 


2.20 


0.375 


0.011 


0.382 


-0.007 


-14.67 


-105.56 


1.32 


51351.2418 


6/22/99 05:48 


2.20 


0.545 


0.023 


0.545 


0.000 


-5.36 


-106.13 


1.84 


51354.2142 


6/25/99 05:08 


2.19 


0.728 


0.028 


0.727 


0.000 


-13.86 


-105.64 


1.37 


51354.2393 


6/25/99 05:44 


2.19 


0.836 


0.038 


0.866 


-0.030 


-2.86 


-106.20 


1.97 


51355.2126 


6/26/99 05:06 


2.18 


0.305 


0.015 


0.291 


0.014 


-13.36 


-105.68 


1.40 


51356.1764 


6/27/99 04:14 


2.20 


0.140 


0.010 


0.147 


-0.007 


-27.47 


-103.85 


0.59 


51356.2139 


6/27/99 05:08 


2.20 


0.240 


0.010 


0.239 


0.001 


-11.59 


-105.81 


1.49 


51356.2534 


6/27/99 06:04 


2.20 


0.375 


0.010 


0.374 


0.001 


5.79 


-106.12 


2.44 


51357.1790 


6/28/99 04:17 


2.20 


0.435 


0.010 


0.439 


-0.003 


-25.31 


-104.22 


0.72 


51357.1932 


6/28/99 04:38 


2.20 


0.417 


0.019 


0.430 


-0.013 


-19.33 


-105.07 


1.06 


51357.2077 


6/28/99 04:59 


2.20 


0.417 


0.011 


0.422 


-0.005 


-13.09 


-105.70 


1.41 
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ABSTRACT 

We report on the determination of the visual orbit of the double-hned spectroscopic 
binary system 12 Bootes with data obtained by the Palomar Testbed Interferometer 
in 1998 and 1999. 12 Boo is a nearly equal-mass double-lined binary system whose 
spectroscopic orbit is well known. We have estimated the visual orbit of 12 Boo from 
our interferometric visibility data fit both separately and in conjunction with archival 
and CORAVEL radial velocity data. Our 12 Boo orbit is in good agreement with 
the spectroscopic results, and the physical parameters implied by a combined fit to 
our visibility data and radial velocity data result in precise component masses. In 
particular, the orbital parallax of the system is determined to be 27.09 it 0.41 mas, 
and masses of the two components are determined to be 1.435 it 0.023 Mq and 1.409 
± 0.020 M0, respectively. 

Somewhat remarkably, even though the two components are nearly equal mass, the 
system exhibits a significant brightness difference between the components in the near 
infrared and visible. We attribute this brightness difference to evolutionary differences 
between the two components in their transition between main sequence and giant 
evolutionary phases, and based on theoretical isochrones we can estimate a system 
age. Further, because the atmospheres of the two components are becoming more 
convective, we suggest the system components are currently at or near synchronous 
rotation, and the system orbit is in the process of circularizing. 



1. Introduction 

12 Bootes (d Bootes, HR 5304, HD 123999) is a sliort-period (9.6 d) binary system with 
nearly-equal mass (q ~ 0.98) components. The system was first detected as a radial velocity 
variable by Campbell & Wright (1900), and the first "good" double- lined orbit was calculated 
by Abt & Levy (1976, hereafter AL76). Merrill (1922) attempted to resolve the system with 
the Mount Wilson interferometer, but did not see any visibility variations and placed it in 
his "Apparently Single" category. Very recently the AL76 orbit has been reconfirmed by an 



independent CORAVEL radial velocity orbit ( pe Medeiros fc Udry 1999| , hereafter DU99). The 



composite system has been consistently assigned the spectral type F8IV - F9IVw, the latter 
by Barry (1970), with the "w" indicating weak ultraviolet metallic features. There is general 
consensus that the components of 12 Boo have evolved off the main sequence. All studies 
seem to confirm that 12 Boo has heavy element abundances near solar proportions ( puncan 
198lt [Balachandran 199q , |Lebre et al.T999| ). 



12 Boo is listed as a triple system by Tokovinin (1997), presumably because the WDS lists 
12 Boo as having a visual companion at a separation of approximately 1" at a position angle of 
8°, but lists no magnitude ( Worley &: Douglass 199^ ). However, McAlister, Hartkopf, & Mason 



(1992) find no companion to 12 Boo within the limits of their speckle observations (separation 
greater than 0.03" and A m < 1.5 mag), and therefore list it as "single" in their Table 5. 

Herein we report the determination of the 12 Boo visual orbit from near-infrared, long-baseline 
interferometric measurements taken with the Palomar Testbed Interferometer (PTI). PTI is a 
110-m H (1.6/im) and i^-band (2.2^m) interferometer located at Palomar Observatory, and 
described in detail elsewhere ( |Colavita et al. 1999| ). PTI has a minimum fringe spacing of roughly 
4 milliarcseconds (10^^ arcseconds, mas) in i^-band at the sky position of 12 Boo, allowing 
resolution of this binary system. We further add photometric and spectroscopic measurements 
in an attempt to understand the fundamental stellar parameters and evolution of the 12 Boo 
components. 



2. Observations 

The interferometric observable used for these measurements is the fringe contrast or visibility 

(squared) of an observed brightness distribution on the sky. Normahzed in the interval [0:1], a 

single star exhibits monochromatic visibility modulus in a uniform disk model given by: 

2J,{nB9/X) 

TrBd/X ^ ' 

where Ji is the first-order Bessel function, B is the projected baseline vector magnitude at the 

star position, 9 is the apparent angular diameter of the star, and A is the wavelength of the 

interferometric observation. The expected squared visibility in a narrow bandpass for a binary 

star such as 12 Boo is given by: 



V^ + Vi r^ + 2ViV2r cos(^ B • s) 



v^, = ^^—^ ^^^_^_~:^^^^^^ ^2) 



where V\ and V2 are the visibility moduli for the two stars separately as given by Eq. Q, r is 
the apparent brightness ratio between the primary and companion, B is the projected baseline 
vector at the system sky position, and s is the primary-secondary angular separation vector on 
the plane of the sky (|Hummel et al. 19951) . The V"^ observables used in our 12 Boo study are 
both narrow-band V'^ from individual spectral channels (Colavita et al. 1999), and a synthetic 
wide-band y^, given by an incoherent SNR- weighted average V"^ of the narrow-band channels in 
the PTI spectrometer ( Colavita 1999 ). In this model the expected wide-band 1^^ observable is 



approximately given by an average of the narrow-band formula over the finite bandpass of the 

spectrometer: 

1 " 
^-^ = -EKV.(A.) (3) 

i 

where the sum runs over the channels covering the infrared //-band (1.5 - 1.8 /im) and i^-band 
(2 - 2.4 /im) of the PTI spectrometer; PTI operating wavebands are excellent matches to the 
CIT photometric system (Colavita et al. 199£, Elias et al. 1982, Elias et al. 1983). Separate 



calibrations and binary model fits to the narrow-band and synthetic wide-band V'^ datasets 
yield statistically consistent results, with the synthetic wide-band data exhibiting superior fit 
performance. Consequently we will present only the results from the synthetic wide-band data. 



Object 


Spectral 


Star 


12 Boo 


Adopted Model 


Name 


Type 


Magnitude 


Separation 


Diameter (mas) 


HD 121107 


G5III 


6.1 V/4.1 K 


8.2° 


0.84 ± 0.06 


HD 128167 


F2 V 


5.9 V/3.9 K 


7.1° 


0.77 ± 0.04 


HD 123612 


K5III 


5.7 V/3.5 K 


0.92° 


1.33 ± 0.10 



Table 1: PTI 12 Boo Calibration Objects Considered in our Analysis. The relevant parameters for 
our three calibration objects are summarized. The apparent diameter values are determined from 
effective temperature and bolometric flux estimates based on archival broad-band photometry, and 
visibility measurements with PTI. 



12 Boo was observed in conjunction with objects in our calibrator list by PTI in i^-band 
(A ~ 2.2/im) on 17 nights between 21 June 1998 and 28 June 1999, covering roughly 39 periods of 
the system. Additionally, 12 Boo was observed by PTI in H-hand (A ~ 1.6/im) on 28 May 1999. 
12 Boo, along with calibration objects, was observed multiple times during each of these nights, 
and each observation, or scan, was approximately 130 sec long. For each scan we computed a 
mean V^ value from the scan data, and the error in the V"^ estimate from the rms internal scatter 
( Colavita 199S| ). 12 Boo was always observed in combination with one or more calibration sources 
within r^ 10° on the sky. For our study we have used three stars as calibration objects: HD 
121107 (G5 HI), HD 128167 (F2 V), and HD 123612 (K5 HI). Table | lists the relevant physical 
parameters for the calibration objects. 

The calibration of 12 Boo V"^ data is performed by estimating the interferometer system 
visibility (V^^^) using calibration sources with model angular diameters, and then normalizing 
the raw 12 Boo visibility by V^y^ to estimate the V'^ measured by an ideal interferometer at that 



epoch ( Mozurkewich et al. 1991 , Boden et al. 199q ). Calibrating our 12 Boo dataset with respect 
to the three calibration objects listed in Table |l] results in a total of 72 calibrated scans (62 in 
K, 10 in H) on 12 Boo over 18 nights in 1998 and 1999. Our calibrated synthetic wide-band V'^ 



measurements are summarized in Table |2|. 

To our PTI visibilities and 17 double-lined AL76 radial velocity measurements we have 
added seven double-lined radial velocity measurements from CORAVEL. These seven CORAVEL 
RV measurements are a subset of the 12 measurements from DU99. The 24 radial velocity 
measurements used in our analysis are summarized in Table ^. 

Finally, on 29 and 30 March 1999 we obtained broad-band infrared photometry on the 12 Boo 
system using the 200" Hale telescope at Palomar Observatory. 12 Boo was observed during 
photometric conditions with respect to Elias IR standards 9529 and 9539 (HD 105601 and HD 
129653 respectively, Elias et al. 1982| ) at airmasses less than 1.2. Table ^ summarizes the results 



of our photometric measurements on 12 Boo. 



3. Orbit Determination 



As in previous papers in this series (Boden et al. 1999a, Boden et al. 1999b| ) the estimation 



of the 12 Boo visual orbit is made by fitting a Keplerian orbit model with visibilities predicted 
by Eqs. |2| and y directly to the calibrated (narrow-band and synthetic wide-band) V'^ data on 



12 Boo (see also Armstrong et al. 1992, Hummel et al. 1993, 1995, 1998). The fit is non-linear in 



the Keplerian orbital elements, and is therefore performed by non-linear least-squares methods 
(i.e. the Marquardt-Levenberg method. Press et al. 1992| ). As such, this fitting procedure takes an 



initial estimate of the orbital elements and other parameters (e.g. component angular diameters, 
brightness ratio), and evolves that model into a new parameter set that best fits the data. 
However, the chi-squared surface has many local minima in addition to the global minimum 
corresponding to the true orbit. Because the Marquardt-Levenberg method strictly follows a 
downhill path in the x^ manifold, it is necessary to thoroughly survey the space of possible binary 
parameters to distinguish between local minima and the true global minimum. In addition, as 
the V'^ observable for the binary (Eqs. and H) is invariant under a rotation of 180°, we cannot 



7 



differentiate between an apparent primary/secondary relative orientation and its mirror image on 
the sky. Consequently there remains a 180° ambiguity in our determination of the longitude of the 
ascending node, 0,, which we quote by convention in the interval [0:180). By similar arguments our 
V'^ observable does not distinguish the longitude of periastron {to) for the primary and secondary 
component. We have constrained our estimate to be grossly (within 180°) consistent with the (Ji 
value of approximately 290° (AL76, DU99). 

In addition to our PTI visibility data we have used the double-lined radial velocity data 
from AL76 and CORAVEL. We incorporate these data into the orbit estimation, utilizing both 



interferometric visibility and radial velocity data either separately or simultaneously ( Hummel 



et al. 199^ , [Boden et al. 1999b| ). The w-degeneracy discussed above is resolved by the inclusion 



of radial velocity data in our orbital solution (Table |5|), however the determination of il remains 
ambiguous by 180°. 

In the case of 12 Boo the parameter space is significantly narrowed by the high-quality 
spectroscopic orbits from AL76 and DU99, and the Hipparcos distance determination sets the 
rough scale of the semi-major axis ([ESA 19*97, Ferryman et al. 1997). Further, at the distance 



of 12 Boo the apparent diameters of the two components of the system are not strongly resolved 
by FTI, so we have constrained the estimated diameters of both components to model values of 
0.63 =b 0.06 and 0.46 ± 0.05 mas for the primary and secondary components respectively (see 
the discussion in § §). Given this limited parameter space, the correct orbit solution is readily 
obtained by exhaustive search for the global minimum in the x^ manifold. 

Figure [l| depicts the relative visual orbit of the 12 Boo system, with the primary component 
rendered at the origin, and the secondary component rendered at periastron. We have indicated 
the phase coverage of our V'^ data on the relative orbit with heavy lines; our data samples most 
phases of the orbit well, leading to a reliable orbit determination. 

Table y lists the complete set of V"^ measurements in our 12 Boo dataset and the prediction 
based on the best-fit orbit model (our "Full-Fit" model. Table g) for 12 Boo. Table ^ gives V'^ 



o 
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Retrograde Orbit 



P = 9.60 d 
a = 3.39 ±0.05 mas 
e = 0.188 ±0.002 
i= 108.6 ± 0.4 deg 
7roj.b = 27.1 ±0.4 mas 




12 Boo Orbit Trace 

V" Phase Coverage 

12 Boo Primary 

12 Boo Secondary at Periastron 




Relative RA (mas) 



Fig. 1. — Visual Orbit of 12 Boo. The relative visual orbit model of 12 Boo is shown, with the 
primary and secondary objects rendered at Tq (periastron). The heavy lines along the relative orbit 
indicate areas where we have orbital phase coverage in our PTl data (they are not separation vector 
estimates); our data sample most phases of the orbit well, leading to a reliable orbit determination. 
Component diameter values are estimated (see discussion in § 0), and are rendered to scale. 



Band 


Magnitude 


^CIT 
KciT 


3.810 ± 0.020 
3.600 ± 0.020 
3.550 ± 0.016 


V (archival) 


4.83 ± 0.01 



Table 4: 12 Boo Near-Infrared Photometric Measurements. We summarize our near-infrared 
photometric measurements on the 12 Boo system (taken at the 200" telescope at Palomar 
Observatory on 29 and 30 March 1999), and archival y-band photometry from the Simbad database. 
Our infrared photometry is taken in the CIT system (Elias et al. 1982, Elias et al. 1983| ). 



measurements and times, measurement errors, model predictions, the photon-weighted average 
wavelength, u — v coordinates, and on-target hour angle for each of our calibrated 12 Boo 
observations. Figures ^ and ^ illustrate our model fit for 12 Boo. Figure ^a shows four consecutive 
nights of PTI V^ data on 12 Boo (28 Feb - 3 Mar 1999), and V^ predictions based on the best-fit 
model for the system (our "Full-Fit" model. Table |5|). Figure ^3 gives a phase plot of V"^ residuals, 
with an inset V'^ residual histogram. The model predictions are in good agreement with the 
observed data, with an rms V'^ residual of 0.033 (average absolute V'^ residual of 0.023), and a 
X^ per Degree of Freedom (DOF) of 1.04. The quality of the V'^ fit is similar to those seen in 
other PTI orbital analyses, and we can see no signs of either bias or excess noise contributed by a 
putative optical companion at a separation of 1" . Figure pi gives a radial velocity phase plot of the 
AL76 and CORAVEL radial velocity data and the predictions of our "Full-Fit" orbital solution. 
The RV rms residual in our solution is 2.3 km s~^ (average absolute residual of 1.7 km s~^), but 
the fit quality of the CORAVEL RV subset (Table 0) is considerably better, with an rms residual 
of 0.59 km s~^ (average absolute residual of 0.44 km s~^). 

Spectroscopic orbit parameters (from AL76 and DU99) and our visual and spectroscopic orbit 
parameters of the 12 Boo system are summarized in Table |5|. We give the results of separate fits 
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Fig. 2. — V"^ Fit of 12 Boo. a: Four consecutive nights (28 February - 3 March 1999) of cahbrated 
y^ data on 12 Boo, and V"^ predictions from the best-fit model for the system. The lower frame 
shows individual V'^ residuals between the calibrated data and best-fit model, b: A phase plot of 
i^-band V"^ fit residuals from our Full-Fit solution (Table ||). We have inset a V'^ error residual 
histogram. 
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Fig. 3.— RV Fit of 12 Boo. A phase plot of radial velocity data from AL76 and CORAVEL and fit 
predictions from our Full-Fit solution (Table |5|). We have inset RV error residual histograms for all 
the RV data, and separately for the DU99 CORAVEL RV; the higher precision of the CORAVEL 
data is evident. 
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to only our V'^ data (our "F^-only Fit" solution), and a simultaneous fit to our V'^ data and the 
double-lined radial velocities from AL76 and CORAVEL (our "Full-Fit" solution) - both with 
component diameters constrained as noted above. For the orbit parameters we have estimated 
from our visibility data we list a total one-sigma error in the parameter estimate, and the 
separate one-sigma errors in the parameter estimates from statistical (measurement uncertainty) 
and systematic error sources. In our analysis the dominant forms of systematic error are: (1) 
uncertainties in the calibrator angular diameters (Table ||); (2) uncertainty in the center-band 
operating wavelength (Aq ~ 2.2 /um), taken to be 20 nm (~1%); (3) the geometrical uncertainty 
in our interferometric baseline ( < 0.01%); and (4) uncertainties in orbital parameters constrained 
in our fitting procedure (i.e. the angular diameters in both solutions, the period in the 'V^-only" 
solution). 

4. Physical Parameters 

Physical parameters derived from our 12 Boo "Full-Fit" visual/spectroscopic orbit are 
summarized in Table ^. As in Table ^, for physical parameters we have estimated we quote total 
one-sigma errors, and statistical and systematic contributions. (Exceptions to this are quantities 
we have estimated using the T^-band spectroscopy discussed in § g the error is taken as statistical, 
and is relatively large compared to the interferometric determinations.) Notable among these 
is the high-precision determination of the component masses for the system, a virtue of the 
precision of the AL76 and CORAVEL radial velocity measurements on both components and the 
moderately high inclination of the orbit. We estimate the masses of the primary and secondary 
components as 1.435 ± 0.023 and 1.408 ± 0.020 Mq, respectively 



The Hipparcos catalog hsts the parallax of 12 Boo as 27.27 ± 0.78 mas ( |ESA 1997|) . The 
distance determination to 12 Boo based on our orbital solution is 36.93 it 0.56 pc, corresponding 
to an orbital parallax of 27.08 it 0.41 mas, consistent with the Hipparcos result at 0.7% and 
0.2-sigma. 
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Orbital 


AL76 


DU99 


PTI 


98/99 


Parameter 


y2.only Fit 


Full Fit 


Period (d) 


9.604538 


9.6046 


5.e{?^5&5 


9.604565 




± 2.2 X 10"^ 


± 1 X 10"'' 




± 1.0 (1.0/0.07) X 10-5 


To (MJD) 


17679.511 


48990.29 


51237.749 


51237.779 




± 0.084 


±0.03 


± 0.024 (0.016/0.018) 


± 0.020 (0.018/0.008) 


e 


0.1933 


0.193 


0.1781 


0.1884 




± 0.0080 


± 0.004 


± 4.8 (4.4/2.0) X 10-3 


± 2.2 (2.2/0.2) X 10-3 


Ki (km s-i) 


67.4 ± 0.8 


67.11 ± 0.41 




67.84 ± 0.31 (0.31/0.03) 


K2 (km s-i) 


66.5 ± 0.9 


70.02 ± 0.48 




69.12 ± 0.48 (0.48/0.03) 


7 (km s^"*") 


9.2 ± 0.4 


9.29 ± 0.19 




9.11 ± 0.13 (0.13/0.01) 


ui (deg) 


290 ±3 


286.19 ± 1.31 


287.0 ± 1.3 (0.9/0.9) 


287.03 ± 0.75 (0.69/0.30) 


^1 (deg) 






9.56 ± 0.41 (0.40/0.07) 


10.17 ± 0.45 (0.40/0.21) 


i (deg) 






108.84 ± 0.32 (0.19/0.26) 


108.58 ± 0.36 (0.29/0.21) 


a (mas) 






3.413 ± 0.039 (0.028/0.027) 


3.392 ± 0.050 (0.036/0.034) 


A/CciT (mag) 






0.614 ± 0.015 (0.011/0.011) 


0.618 ± 0.022 (0.019/0.011) 


Aii'ciT (mag) 






0.588 ± 0.066 (0.063/0.019) 


0.566 ± 0.066 (0.063/0.019) 


x7dof 


1.4 


0.90 


0.82 
0.023 


1.2 (1.04 1/2/2.8 RV) 


\Rv^ 


0.023 


Rrv\ (kms-1) 


1.7 (2.2 AL/0.44 COR) 



Table 5: Orbital Parameters for 12 Boo. Summarized here are the apparent orbital parameters for 
the 12 Boo system as determined by AL76, DU99, and PTI. We give two separate fits to our data, 
with and without including the double-lined AL76 and CORAVEL radial velocities in the fit. For 
parameters we have estimated by including our PTI observations we separately quote one-sigma 
errors from both statistical and systematic sources, given as {(Jstat/<^sys)-, and the total error as the 
sum of the two in quadrature. Quantities given in italics are constrained to the listed values in 
our model fits. We have quoted the longitude of the ascending node parameter (fi) as the angle 
between local East and the orbital line of nodes measured positive in the direction of local North. 
Due to the degeneracy in our V"^ observable there is a 180° ambiguity in i7; by convention we quote 



it in the interval of [0:180). We quote mean absolute V and RV residuals in the fits, |-Ry2| and 
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Component Diameters and Effective Temperatures The "effective" net angular diameter 
of tlie 12 Boo system has been estimated using the infrared flux method (IRFM) by Blackwell 
and collaborators ( [Blackwell et al. 199C , Blackwell Sz Lynas-Gray 1994| ) at approximately 0.8 



mas. At this size neither of the 12 Boo components are resolved by PTI, and we must resort to 
model diameters for the components. Following Blackwell, we have estimated 12 Boo component 
diameters through bolometric flux and effective temperature (Tefj) arguments. Blackwell and 
Lynas-Gray (1994) list the bolometric flux of the 12 Boo system at 3.11x10^^'^ W m^^, and 
Teff as 6204 K, both without error estimates. By assuming our A'-band flux ratio as a surrogate 
for the bolometric flux ratio (correct in the limit of similar T^s for the two components) it is 
straightforward to compute component diameters as a function of individual component Tefj. If we 
assume the 6204 K Tgff for both components we arrive at diameter estimates of 0.62 and 0.46 mas 
for the primary and secondary components respectively. However, comparing our interferometric 
component magnitude difference in i^-band and a spectroscopic estimate in V-band as we shall 
see in § 1^ we can employ an empirical effective temperature - {V - Kn) color index relationship 
for subgiant stars derived by Blackwell et al. (1990) to estimate Teg of 6050 it 75 and 6250 it 
130 K for the more evolved primary and less-evolved secondary components respectively. (In this 



computation we have corrected for the increase in the 12 Boo Tpg estimate between Blackwell et 



al. 1990| and Blackwell &: Lynas-Gray 1994, and have used the {V - K) to {V - Kn) empirical 



color index correction of Selby et al. 198S with the component {V - K) values of Table g.) With 
these individual component Tgfj we find model diameters of 0.65 and 0.46 mas for the primary 
and secondary components respectively, similar to the values derived from applying the Blackwell 
Teff to both components. We have averaged these two sets of diameter estimates for our final 
model diameters, arriving at 0.63 and 0.46 mas for the primary and secondary component angular 
diameters respectively. Rigorous error estimates in these diameter models are made impossible 
without estimates of the errors in the input quantities; we have taken ad hoc 10% one-sigma errors 
in these models in our orbit systematic error calculations. At the consensus distance estimate to 
12 Boo these model angular diameters correspond to model component linear radii of 2.51 ± 0.25 
and 1.83 ± 0.18 Rq for the primary and secondary components respectively. These linear radii 
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values are roughly a factor of two smaller than the putative Roche lobe radii for these two stars 



(Iben 1991 Eq. 1), making significant mass transfer unlikely at this stage of system evolution. We 
have further ignored any corrections due to stellar limb darkening; at these putative component 
sizes our data would be highly insensitive to limb-darkening effects (for examples of observational 
consequences of stellar limb darkening see Quirrenbach et al. 1996| , pajian et al. 1998 , and 
references therein). 

5. Component Intensity Ratio 

As given in Table ^, despite the fact that the mass ratio of 12 Boo is near 1 (q ~ 0.98), our 
best fit model of the system requires that the near-infrared luminosities of the two components 
are significantly different. This brightness difference appears to be verified in the visible in 
high-resolution spectroscopy (see ^. 

The most sensitive test of this relative component brightness difference is the V'^ value at 
interferometric visibility minima; following Eq. |2|, for unresolved components the minimum V'^ on 
a binary system is approximated by: 

v^- ^('-^y (4) 

In the limit that r -^ 1, y^ ~ 0, and a fringe tracking interferometer like PTI cannot track the 
source ( Colavita et al. 199S| ). Having seen indications of the intensity asymmetry in our initial 



orbit solution, on 15 April 1999 we performed a deliberate experiment to measure the V'^ on 
12 Boo through a predicted visibility minimum. The results are given in Figure §. Both raw and 
calibrated /C-band V'^ measurements on 12 Boo are shown, along with the predictions derived 
from our "Full-Fit" orbit model with a 57 it 1 % component intensity ratio (Table ^), and the 
same orbit geometry but an assumed 90% component intensity ratio. The fact that the raw and 
calibrated V'^ measurements are significantly above the 90% intensity ratio model is unequivocal 
evidence that the K-hand component brightness ratio is significantly less that unity. Moreover, 
our best-fit orbit model in general, and our X-band intensity ratio estimate (0.57, Table y) in 
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Physical 
Parameter 



Primary 
Component 



Secondary 
Component 



a (10-2 AU) 

Mass (Mq) 

Sp Type (Barry 1970) 

System Distance (pc) 

TTorb (mas) 

Model Diameter (mas) 

Mi^_ciT (mag) 

Mh-cit (mag) 

My (mag) 

V-K (mag) 



6.205 ± 0.032 (0.031/0.008) 6.322 ± 0.046 (0.045/0.008) 
1.435 ± 0.023 (0.023/0.004) 1.408 ± 0.020 (0.019/0.004) 

F9 IVw 

36.93 ± 0.56 (0.43/0.37) 

27.08 ± 0.41 (0.31/0.27) 

0.63 {± 0.06) 0.46 {± 0.05) 

1.200 ± 0.038 (0.031/0.027) 1.818 ± 0.039 (0.032/0.029) 

1.269 ± 0.048 (0.037/0.034) 1.835 ± 0.063 (0.048/0.044) 

2.524 ± 0.052 3.024 ± 0.077 

1.324 ± 0.044 1.206 ± 0.072 



Table 6: Physical Parameters for 12 Boo. Summarized here are the physical parameters for the 
12 Boo system as derived primarily from the Full-Fit solution orbital parameters in Table |^. 
Quantities listed in italics (i.e. the component diameters, see text discussion) are constrained to the 
listed values in our model fits. As for all our PTI-derived orbital parameters we have quoted both 
total error and separate contributions from statistical and systematic sources (given as o-gtat/f^sys), 
with the exceptions of quantities involving the spectroscopic AV determination where the error 
is assumed to be statistical, and is large compared with the interferometric estimates. Infrared 
absolute magnitudes are quoted in the CIT system ([Elias et al. 198^, Elias et al. 1983). 
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particular, does an excellent job of predicting the V"^ variations (Eq. ^) and the minimum V"^ 
(Eq. ^ respectively. 

In 1999 PTI added the capability to make if -band (A ~ 1.6 /um) visibility measurements. As 
mentioned in § ^, 12 Boo was observed with PTI in //-band on one night (28 May 1999) from 
which we estimate an ii^-band magnitude difference of the two components. Figure ^ shows the 
calibrated H-h&nd visibilities obtained on 12 Boo, and a priori (using the i^-band intensity ratio) 
and fit predictions for 12 Boo. This small amount of H-h&nd data indicates that the magnitude 
difference of the 12 Boo components in H is 0.566 =b 0.045 (Table ^); the difference between this 
and our i^T-band value of 0.618 ± 0.021 is not formally significant. 

Comparisons With Stellar Models Given our estimates of component masses, absolute 
magnitudes, and color indices derived from our measurements and orbital solution (Table y), we 
can examine the 12 Boo components in the context of stellar models. Figure |^ shows the placement 
of the 12 Boo components in observable parameter mass-magnitude and color-magnitude spaces, 
along with theoretical isochrone evolutionary tracks from Bertelli et al. (1994, hereafter B94). 
Here we have used B94 isochrone tracks for Z = 0.02 (B94 Table 5) based on the assumption 



of solar abundances for the two 12 Boo components. The [Fe/H] values given in ( Duncan 



198H iBalachandran 199C| , [Lebre et al. 1999|) indicate Z12 boo = 0.0177 ± 0.0017 (B94 Eq. 10), 



hence the applicability of the solar abundance models. The isochrones suggest that the primary 
component of 12 Boo is in the midst of particularly rapid evolution; unfortunately the coverage of 
the isochrone models is relatively coarse in the region of interest for this system. B94 uses Johnson 
infrared passbands in their calculations; for comparison with B94 models we have transformed 
our component Adx values (Table 0) from the CIT to the Johnson system using the relation from 
Bessell & Brett (1988). 

By jointly fitting our component My, Mk, and V - K estimates (Table ^ to the B94 model 
predictions of these quantities at our estimated component masses as a function of age we find a 
best-fit 12 Boo system age estimate of 2.550 ± 0.023 Gyrs. The quoted error in this age estimate is 
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Fig. 4. — 12 Boo Visibility Scan Through Minimum. In order to test the brightness asymmetry 
hypothesis, a tentative 12 Boo orbit model was used to predict a V"^ minimum on 15 April 1999. 
Both raw and calibrated if -band V"^ measurements on 12 Boo are shown, along with the predictions 
derived from our "Full-Fit" orbit model with a 57% component intensity ratio (Table ^), and the 
same orbit but an assumed 90% component intensity ratio. The fact that the raw and calibrated 
V'^ measurements are significantly above the 90% intensity ratio model is unequivocal evidence 
that the K-hand component brightness ratio is significantly less that unity. 
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Fig. 5. — 12 Boo H-hand Visibilities. In order to estimate the differential i:f -magnitude of the two 
12 Boo components we observed 12 Boo on one of our first H-hand science commissioning runs at 
PTl on 28 May 1999. Our i^-band derived 12 Boo orbit model, including the iT-band component 
intensity ratio, was used to predict the expected V"^ variations at f/^-band. A fit to the H-hand data 
for component intensity ratio results in a slightly smaller component magnitude difference than we 
derive in our iC-band fit, but the difference is not formally significant (see text and Table 0). 
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purely statistical and does not include possible systematic error effects in the B94 isochrones. The 
fit is of relatively poor quality, with a chi-squared of 10.3 for two degrees of freedom. If instead we 
use only the observed component magnitude differences in V and K (Table y) compared to the B94 
models we find a best-fit system age of 3.220 it 0.022 Gyrs (again the error is purely statistical). 
That these two different age estimates do not agree within their statistical errors indicates a 
significant discrepancy between the B94 models and our observed parameters for the components 
in the 12 Boo system. This discrepancy makes it difficult to decide which of the two different 
estimates is more accurate. A simple average of these two different age estimates would suggest a 
12 Boo system age of 2.89 it 0.36 Gyrs; the error in this age estimate is apparently dominated by 
systematic errors in the application of the B94 stellar models to the 12 Boo components. 

6. Tidal Interaction and Component Rotation in 12 Boo 

In short-period binary systems, the components gravitationally interact so as to circularize 



the orbit and synchronize the component rotations to the orbit period ( Zahn 1977 , Hut 1981 ). 
In a survey of nearby "solar-like stars", Duquennoy & Mayor (1991) give evidence that systems 
with periods shorter than 11 days are "all circularized due to tidal effects occurring during their 
evolution on the main sequence." The circularization and synchronization phenomena necessarily 
require an energy dissipation mechanism, which is generally thought to be associated with 



convection in the outer envelopes of evolved stars (Verbunt &; Phinney 1995). 



12 Boo is interesting from a tidal interaction perspective because the system orbit is modestly 
eccentric (e ~ 0.19; Table P), yet the components of the system have in fact evolved off the main 



sequence. With our masses (Table p[) and eclipsing binary-derived models (Popper 1980, Andersen 



1991) we infer that during their main-sequence lives the 12 Boo components were of approximate 
spectral type F3 V. This typing would suggest the reason for the remnant eccentricity in the 
12 Boo system is the putative lack of a convective dissipation mechanism in early-F main sequence 
stars. However, as the components evolve their atmospheres become more convective, and tidal 
circularization and synchronization should begin. 
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Fig. 6. — 12 Boo Components In Mass-Magnitude and Color-Magnitude Spaces. Here we depict the 
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While the orbit remains eccentric, the characteristic timescale for rotational synchronization 
in close binaries with convective envelopes is approximately given by (Zahn 1977], Giuricin et 



al. 1984D : 

fR\'^ -2 

4ynch ~ ( ~ ) 9 yr 

with R as the stellar radius, a as the orbital semi-major axis, and q as the binary mass fraction. 
Our models for the 12 Boo components and orbit imply synchronization timescales of 1.6x10^ 
and 1.1x10^ yr for the primary and secondary components respectively. Because these are 
significantly shorter than the likely system age of ~ 3 Gyr (§ ^) it is interesting to look for signs 
of synchronization in 12 Boo. 



Several recent measurements of the rotation vsini of 12 Boo exist (Balachandran 199C, De 
Medeiros et al. 1997 , [Lebre et al. 1999| , DU99), offering the possibility to test whether the two 
components are synchronously rotating. Additionally we have taken a high resolution ELODIE 
spectrum of the system to assess the rotation and the relative brightness of the components (see 
below). We summarize these recent rotation measurements and our own in Table 0. The consensus 
is that both 12 Boo components are rotating considerably faster than the mean rotational velocity 
of 5.4 km s~^ for subgiant stars with (B-V) ~ 0.55 (De Medeiros et al. 1996|). 



ELODIE Measurements In June 1999 we took a spectrum of the 12 Boo system with the 
ELODIE high resolution echelle spectrograph at the Observatoire de Haute Provence, in France 
( Baranne et al. 1996| ). A cross-correlation function (CCF) has been computed using the on-line 



reduction package available with this instrument (Figure |^ . 

The CCF corresponds to a mean of the spectral lines selected by the template ( iQueloz 1994 ) 



The width of the CCF is a measurement of the mean broadening of spectral lines. It can be 
calibrated for each spectral type to provide a vsini measurement of the star ( Queloz et al. 19*9^ ). 



The equivalent width of the CCF (W) is sensitive to the temperature and the metallicity of the 
star. If the stellar spectral lines have no strong asymmetries or the star is not a fast rotator the 
shape of the CCF can be very well approximated by a Gaussian function {G{v)) in absorption. 
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Primary 


Secondary 




vsini 


vsini 




(km s~i) 


(km s-i) 




10 ± 3 

12.7 ± 1 

12.7 ± 1 

12.5 (± 1) 




Palachandran 1990 








De Medeiros et al. 1997 








Lebre et al. 1999| 




DU99 


9.5 (± 1) 


ELODIE/This Work 


13.1 ± 0.3 


10.4 ± 0.3 



Composite 



13.00 ± 0.28 10.33 ± 0.29 



Model Synch Rotation 12.5 (± 1.2) 9.1 (± 0.9) 

Model Pseudo-Synch Rotation 15.2 (± 1.5) 11.1 (± 1.1) 

Table 7: Recent vsini Measurements for 12 Boo Components. Summarized here are the most 
recent vsini measurements for the 12 Boo system components. For references where a single 
vsini measurement is listed we have assumed this pertains to the primary component. DU99 does 
not list errors for their component vsini estimates; we have arbitrarily taken 1 km s^^ so as to be 
consistent with the characteristic accuracies of earlier CORAVEL determinations (see discussions 
De Medeiros et al. 1996| , pe Medeiros et al. 1997| , [Lebre et al. 1999|) . We list weighted average 



m 



V sin i estimates based on the listed measurements. For comparison we further give model estimates 
of vsini for synchronous and pseudo-synchronous rotation of the two components assuming the 
physical sizes discussed in § ^; errors are based on ad-hoc estimates of systematic errors in the 
diameter models. 
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Fig. 7. — Cross-correlation function of the 12 Boo ELODIE spectrum (filled dots). We see the 
double line feature characteristic of a double-lined spectroscopic binary system. The superimposed 
dotted line is the double Gaussian fit. Empty dots at the bottom display the residuals (cr = 0.07% 
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(l-Giv)). 

In a double-lined spectroscopic binary such as 12 Boo, the equivalent width ratio 
{W1/W2) can be used as a flux ratio indicator if the intrinsic equivalent width ratio 
((1^1/1^2)™*"'^"^) can be estimated from two similar single spectra: Wi = Wj™*"'^"^Fj(X]j -F'j)"^, 
then W1/W2 = (VFi/W2)™*'"^''^ {F1/F2). The mass of each component in the 12 Boo system is very 
similar (Table ^. The CCF equivalent width ratio computed using respectively only the blue 
and the red part of the ELODIE spectrum show no significant differences, suggesting very similar 
Teff. Since the system is coeval a metallicity difference between the two objects seems unlikely. 
Therefore we make the assumption that the intrinsic equivalent width of each of the spectral 
line systems is similar, and use our computed CCF as an indicator of relative brightness. With 
this assumption a magnitude difference Am^ = 0.5 it 0.1 is found, in reasonable agreement with 
the component magnitude difference seen in the infrared interferometric data (§ 1^. The derived 
intrinsic equivalent width agrees with a {B — V) ^ 0.55 star with solar metallicity. Note that the 
gravity difference expected between the two stars is too small to have a significant impact on the 
spectral line equivalent widths. 

To compute the vsini of each of the stars from the CCF, the Queloz et al. (1998) calibration 
has been used. This leads to vsmi= 13.1 ± 0.3 km s^^ for the bright (primary) component and 
vsmi= 10.4 lb 0.3 km s^^ for the dimmer (secondary) component (Table 1^. Note that the smaller 
gravity of the brighter component compared to the other one desaturates some lines from the 
bright object. This would cause an underestimate of vsini, however the effect is small in our 
analysis, and well within the 0.3 km s~^ error in our measurements. 

Assuming the 12 Boo component model diameters calculated above (§ ^), and coplanarity of 
orbital motion and component spin (required for tidal equilibrium), we can calculate the expected 
vsini speeds for synchronous or pseudo-synchronous rotation (where the rotation equilibrates 
with the orbital motion at periapsis when the tidal forces are at a maximum ( [Hut 1981 )) with 



the orbital motion; these are listed in Table |^. The quoted errors in these rotation calculations 
are dominated by the ad hoc uncertainties in our component model diameters. Our ELODIE 
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V sin i measurements suggest that both components are rotating shghtly faster than synchronously, 
but neither component's vsvui seems consistent with the possibihty of pseudo-synchronous 
rotation. However, we would add that the DU99 vsvui values lead one to conclude that both 
components are in essentially synchronous rotation. Because DU99 do not list errors in their 
vsvui values it is difficult to assess the relation between our measurements and theirs. 

Finally, in an attempt to independently confirm component rotation periods through 
differential photometry we consulted the Automated Astronomy Group at Tennessee State 
University. Both the Phoenix-10 automatic photoelectric telescope (APT) in Phoenix, AZ, and 
the Vanderbilt/ Tennessee State 16-inch APT at Fairborn Observatory near Washington Camp, 
AZ, have routinely observed 12 Boo as a photometric reference star, and thus have a large body 
of differential photometric data on the system. They report F-band night-to-night scatter at 
the 0.004 mag level (consistent with scatter on other photometric references), and no discernible 
modulation associated with component rotation ( [Henry 19991 ). 

7. Summary and Discussion 

By virtue of our interferometric resolution and the precision of the AL76 and CORAVEL 
radial velocity data we are able to determine accurate physical parameters for the 12 Boo 
constituents, and an accurate system distance. Our 12 Boo distance estimate is in excellent 
agreement with the Hipparcos trigonometric determination. Our finding of unexpectedly large 
relative K, H, and ^-magnitude differences in the two nearly-equal mass 12 Boo components 
suggests that the system is in a unique evolutionary state, with both components making the 
transition off the main sequence. 

The agreement between our findings on the physical parameters of the two 12 Boo components 
and the stellar models of B94 is not particularly good. The fundamental reason for the discrepancy 
is the unexpectedly large magnitude difference we see in the interferometric and spectroscopic 
data; clearly it leads one to question the veracity of the data in this regard. As described in § |5|, 
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for the interferometric data we performed dedicated observations looking to verify the component 
brightness asymmetry. We have further looked in detail at the individual spectral channel data 
for both the standard PTI pipeline processing, and for alternate detector calibrations ( |Colavita 



1999). In all cases the conclusion is the same: the interferometric visibility measurements 
remain significantly above zero through the visibility minima. Unless there is some unknown and 
heretofore unseen positive bias in our visibility processing at low visibilities (such effects should 



have been evident in Boden et al. 1999b| in particular), the conclusion of a component brightness 



asymmetry seems sound. Having seen the brightness difference in the near-infrared, we verified 
this conclusion with visible spectroscopy as discussed in § g. 

Accepting the component brightness difference as genuine, and further assuming that the 
component metallicities are similar and the stars coeval leads us to the conclusion that the 
dominant source of error in the age estimate for the 12 Boo system arises in B94's treatment of 



stars at this stage of evolution. The system metallicity is reasonably well-determined (Duncan 



1981, Balachandran 1990| , Lebre et al. 1999| ), and our spectroscopy shows no indications of 



dissimilar metallicities between the components. While one might expect possible metallicity 
differences in an interacting binary system, computations indicate that neither component of the 
12 Boo system has filled its Roche lobe (§ ^, making it unlikely that significant mass transfer has 
occurred to date in the system. The calculated separation of the two components is approximately 
0.12 AU, a strong suggestion that these stars formed from the same molecular cloud. Finally, 
there is no indication in the IRAS mid-IR photometric measurements for the 12 Boo system of an 
infrared excess, which might indicate the presence of dust (as in a common envelope or mass- loss 
from one of the components). 

We believe the two most likely reasons for the discrepancy between our component parameters 
(Table 0) and the B94 model tracks stem from the B94 treatment of convection. Iben (1991) 
indicates that the time during which a 1.5 Mq star initially moves to the right of the main 
sequence on an H-R diagram (before the hook back to the left) is the time when hydrogen is 
burning in a convective core. Near the approximate derived age of the 12 Boo system (log T = 
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9.4 yr, 2.51 Gyr), B94 note that there is a dip in the relation between the absolute magnitude and 
age for the turn-off from the main sequence (at the point where the star transitions from radiative 
to convective core H-burning). Further, while B94 do treat overshooting from convective cores in 
their models, only two empirical values (based on measurements of galactic clusters) are used for 
Ac, 0.25 for 1.0 Mq < M < 1.5 M© and 0.5 for M > 1.6 Mq. The fact that the 12 Boo system 
age is likely close to log T = 9.4 yrs, and that the component masses are close to the 1.5 Mq 
B94 transition point in Ac is suggestive that the B94 models may be too coarse in this region 
of parameter space to describe the particular stage of component evolution seen in the 12 Boo 
system. Conversely, we are optimistic that the physical parameters derived here for the 12 Boo 
components are of sufficient quality to service as significant constraints to evolutionary modeling. 

One potential additional contributor to the data-model discrepancy is the possibility of 
biases in our mass estimates. While we see good statistical agreement between the DU99 and 
our "Full-Fit" solution (e.g. 1.2-sigma in the component semi-amplitudes), our solution calls for 
a smaller component mass difference than the DU99 solution would indicate. For example, were 
we simply to adopt the DU99 semi-amplitudes and our "l/^-only" astrometric orbit we would 
obtain component mass estimates of 1.464 it 0.031 and 1.403 ± 0.028 Mq for the primary and 
secondary components respectively. These alternative mass values are statistically consistent with 
our "Full-Fit" values (Table 0), and would alleviate some, but not all, of the discrepancy with the 
B94 models (see Figure ^a in particular). The reason for the potential mass-bias in our results is 
clear: our inclusion of the AL76 RV data, which led AL76 to in fact conclude that the 12 Boo 
secondary was more massive than the primary. 

Although the short-period orbit of 12 Boo has not yet circularized, tidal interactions should 
be at work in the system as both components become more convective with their evolution off 
the main sequence. The body of v sin i measurements suggests that both components of 12 Boo 
are rotating at or slightly higher than synchronous rotational speed, and significantly faster than 



the mean rotational speed for similar subgiant stars ( pe Medeiros et al. 19961 ). Taken at face 
value the v sin i measurements suggest that the primary is closer to synchronous rotation than the 
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secondary; this seems consistent with the rough factor of 7 difference in the synchronization time 
scale expected from the model size difference between the two components. 

The work described in this paper was performed at the Infrared Processing and Analysis 
Center, California Institute of Technology, and the Jet Propulsion Laboratory under contract 
with the National Aeronautics and Space Administration. Interferometer data were obtained 
at Palomar Observatory using the NASA Palomar Testbed Interferometer, supported by NASA 
contracts to the Jet Propulsion Laboratory. Science operations with PTI are conducted through 
the efforts of the PTI Collaboration ( http://huey.jpl.nasa.gov/palomar/ptimembers.html ), and we 
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TABLE 3 

12 Boo RV Dataset, "Full-Fit" Model Predictions, Residuals, and Phases. 



MJD 



Prim 
RV 



Prim 
Error 



Sec 
RV 



Sec 
Error 



Prim 
Model 



Sec 
Model 



Prim 
Resid 



Sec 
Resid 



Orbit 
Phase 



Source 



39247.415 


-33.8 


0.5 


46.9 


2.6 


-32.34 


51.34 


-1.46 


-4.44 


0.60 


AL76 


39277.272 


-53.9 


0.5 


62.3 


2.2 


-51.70 


71.06 


-2.21 


-8.76 


0.71 


AL76 


39317.207 


-43.3 


0.9 


56.6 


0.8 


-39.72 


58.86 


-3.58 


-2.26 


0.86 


AL76 


39909.525 


-13.5 


1.1 


32.8 


1.6 


-16.54 


35.24 


3.03 


-2.44 


0.53 


AL76 


39929.414 


-32.0 


0.9 


52.0 


2.7 


-34.12 


53.15 


2.11 


-1.15 


0.61 


AL76 


39984.255 


46.1 


0.8 


-32.0 


2.0 


48.48 


-31.01 


-2.38 


-0.99 


0.32 


AL76 


39987.246 


-39.5 


0.6 


56.1 


1.4 


-38.72 


57.84 


-0.78 


-1.74 


0.63 


AL76 


40281.453 


64.4 


0.4 


-43.7 


1.0 


64.08 


-46.90 


0.32 


3.20 


0.26 


AL76 


40373.202 


-48.1 


0.8 


73.6 


1.0 


-51.80 


71.17 


3.70 


2.43 


0.81 


AL76 


40374.259 


-13.4 


1.1 


35.5 


1.2 


-14.82 


33.49 


1.42 


2.01 


0.92 


AL76 


40375.247 


46.0 


1.4 


-27.1 


2.0 


46.75 


-29.24 


-0.75 


2.14 


0.02 


AL76 


40613.460 


-46.0 


0.8 


71.7 


1.1 


-49.38 


68.71 


3.38 


2.99 


0.83 


AL76 


40615.493 


53.6 


0.8 


-32.1 


1.6 


53.91 


-36.54 


-0.31 


4.44 


0.04 


AL76 


40666.421 


42.6 


1.0 


-24.4 


1.5 


40.98 


-23.36 


1.62 


-1.04 


0.34 


AL76 


40695.345 


37.0 


0.9 


-20.5 


2.0 


37.49 


-19.81 


-0.49 


-0.69 


0.35 


AL76 


40697.342 


-22.8 


0.9 


39.4 


1.9 


-23.16 


41.99 


0.36 


-2.59 


0.56 


AL76 


41025.473 


-52.5 


0.8 


69.4 


1.6 


-53.39 


72.79 


0.89 


-3.39 


0.72 


AL76 


46957.865 


26.48 


0.26 


-8.18 


0.38 


26.55 


-8.65 


-0.07 


0.47 


0.39 


DU99 


48706.103 


18.44 


0.26 


-2.58 


0.39 


19.92 


-1.90 


-1.48 


-0.68 


0.41 


DU99 


48725.112 


26.6 


0.29 


-8.53 


0.49 


26.30 


-8.41 


0.30 


-0.12 


0.39 


DU99 


48748.975 


-36.6 


0.17 


56.32 


0.25 


-36.65 


55.74 


0.05 


0.58 


0.87 


DU99 


48750.058 


23.92 


0.45 


-5.34 


0.65 


23.91 


-5.97 


0.01 


0.63 


0.99 


DU99 


48753.124 


51.67 


0.59 


-34.3 


0.74 


51.50 


-34.08 


0.17 


-0.22 


0.30 


DU99 


48753.975 


25.79 


0.32 


-6.55 


0.55 


24.73 


-6.81 


1.06 


0.26 


0.39 


DU99 



Note. — Units of all RV in km s 



